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ABSTRACT 



In an attempt to understand the source and propagation of galactic cosmic rays we 
have employed the Modified Weighted Slab technique along with recent values of the 
relevant cross sections to compute primary to secondary ratios including B/C and Sub- 
Fe/Fe for different galactic propagation models. The models that we have considered are 
the disk-halo diffusion model, the dynamical halo wind model, the turbulent diffusion 
model and a model with minimal reacceleration. The modified weighted slab technique 
will be briefly discussed and a more detailed description of the models will be given. We 
will also discuss the impact that the various models have on the problem of anisotropy 
at high energy and discuss what properties of a particular model bear on this issue. 
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Subject headings: Cosmic rays: general — magnetic fields — diffusion — particle ac- 
celeration — shock waves 

1. The application of the Modified Weighted Slab technique to simplified models 

The weighted slab technique has long been used in studying the propagation of cosmic rays 
in the galaxy from their points of origin to their observation points near the earth (Davis 1960; 
Ginzburg & Syrovatskii 1964; Ginzburg & Ptuskin 1976; Lezniak 1979), also see Webber (1997, 
and therein). Several approximations are used in deriving this technique among them the assump- 
tion that energy loss and/or gain is not significant and that the propagation in and loss from the 
galaxy may be described by a function of energy per nucleon alone. Both of these simplifications 
are known to be untrue, for low energies ionization energy loss can be significant and rigidity, or 
energy per charge, is believed to be the parameter that best describes propagation. 

Ptuskin, Jones & Ormes (1996) showed how the weighted slab technique could be made exact 
for galactic propagation models in which energy gains and losses were proportional to the same 
mass density that determined nuclear fragmentation and time dependent processes, eg. radioactive 
decay, do not play a role. This modification allows for the fact that particles had different (usually 
higher) energies in the past and hence different propagation properties and that propagation is 
considered to be a function of rigidity although energy per nucleon is the proper parameter for 
nuclear fragmentation calculations. Strictly, this technique is rigorous only for models in which the 
particle propagation parameters are proportional to a single function of energy for each particle 
species, and hence does not apply to galactic wind models or turbulent diffusion models. However, 
most of these models may be closely approximated by simplified homogeneous models in which the 
mean path length has an exponential distribution with a mean path length that is a particular 
function of rigidity. It is models of this type and approximation that we discuss in this study. 
In this paper we present some results of the numerical simulations where the most recent set of 
spallation cross sections was used. 

It should be noted that these models bear a similarity to the well known leaky box model in that 
they all (at least approximately) yield an exponentially decreasing path length distribution. They 
differ in the manner in which the mean path length varies as a function of energy. The similarity ends 
here, however, as the leaky box model can not predict any anisotropy as it considers the cosmic rays 
to be homogeneously distributed in a "box" of indeterminate size and shape. Furthermore, the only 
independent parameter in the leaky box model is the amount of matter traversed or "grammage" 
which is sufficient for stable nuclei but radioactive nuclei decay as a function of time not grammage 
so the more realistic models can predict significantly different results for these unstable particles. 
Since we will be dealing with stable nuclei only in this work our mathematics will look very much 
like calculations with the leaky box model but it should be remembered that we have significantly 
different physical models in mind when we do these calculations. 
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2. Cross sections 

The cross sections used here now include a completely updated cross section file for the propa- 
gation program. This includes the new primary cross sections in hydrogen targets for C through Ni 
at 600 Mev/n as described in Webber, et al. (1998ab), as well as the hydrogen cross sections for 
essentially all of the secondary nuclei from Li through Mn also at 600 Mev/n reported in Webber 
et al. (1998c). The energy dependence of these isotopic cross sections is updated and extended 
as well using our earlier charge changing cross sections measured between 300 and 1700 Mev/n 
(Webber et al. 1990) and at 15 GeV/n (Webber et al. 1994) and assuming that the isotopic 
fractions are generally independent of energy as confirmed by these earlier measurements and those 
of the Transport Collaboration (Chen et al. 1997). 



3. Disk-halo diffusion model. 
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Fig. 1. — Simplified galactic model with a diffusive halo, galactic matter disk is of infinitesimal 
thickness. 



This propagation model, (Ginzburg, Khazan & Ptuskin 1980; Berezinskii et al. 1990)schemat- 
ically shown in Figure (1), involves a thin (approximately infinitesimal) disk of matter and an ex- 
tended, matter free, halo through which the cosmic rays diffuse with a rigidity dependent diffusion 
coefficient, D(R). The distribution function f(z,p) normalized as N = Air f dpp 2 f (where N is the 
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total cosmic ray number density) obeys the equation 



d df fiva 



dz dz 



1 d 



m 



qo(p)S(z) 



(3-1) 



Here D(p, z) is the particle diffusion coefficient, \x ~ 2.4 mg/cm 2 (Ferriere 1998) is the surface 
mass density of the galactic disk, v = /3c is the particle velocity, a is the total spallation cross section, 
m is the mean mass of interstellar atom, (dp/dt) ion = (i5(z)bo(p)/m < describes the ionization 
energy losses and qo(p)5(z) is the source term that may include the yield from the fragmentation 
of more heavy nuclei. 

We shall assume that diffusion does not depend on position, i.e. D = D(p). There is a cosmic 
ray halo boundary at i z i= H where cosmic rays freely exit from the Galaxy. 

Integrating equation (3-1) in the vicinity of galactic plane lim Jt £ dz{...) at e — ► one can 
find the boundary condition at z = + e : 



9n dfo »va 1 d 
oz m p z op 
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P — fo 
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Qo(p), 



(3-2) 



where fo(p) = f(z = 0,p) is the distribution function at the galactic midplane. 

The solution of equation (3-1) at z ^ under the boundary condition /(i z \= H) = is: 

H — i z i 



/ = /o : 



H 



(3-3) 



Calculating —D^fr from equation (3-3) and substituting it in equation (3-2), one can get the 



following closed equation for fo(p) : 



2Df | a | 1 d 
fivH m p 2 dp 



P — fo 

mv 
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(3-4) 



Now we introduce cosmic ray intensity as function of kinetic energy per nucleon I{E} t )dEf : = 
vfo(p)p 2 dp, so that fo(p) = I(Ek)A/p 2 (A is the atomic number). Equation (3-4) gives the following 
equation for I(Ek) : 



/ d 

+ 



dEi 



I 



X di f dE k [\ dx 
with the effective escape length for diffusing particles 
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(3-5) 



(3-6) 



and the source term 
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Q = *M?!i (3-7) 



Equation (3-5) is identical to the transport equation for cosmic rays in the so called leaky box 
model, the popular empirical model used in the studies of cosmic ray propagation. The leaky box 
model has the only parameter, the escape length Xn,, which describes the cosmic ray propagation 
and this parameter is considered as an empirical characteristic determined from cosmic ray data. 
The physical meaning of the escape length is that it is equal to the mean thickness of matter 
traversed by cosmic rays before exit from the Galaxy. 

The foregoing consideration confirms the well known result that the diffusion model with 
relatively thin galactic disk (with half thickness h <C H) and with flat halo (radius of the galactic 
disk R 3> H) is almost equivalent to the leaky box model for calculation of abundance of stable 
nuclei in cosmic ray (Berezinskii et al. 1990; Ptuskin et al. 1997). The equivalence holds for not 
very heavy nuclei which have the total cross sections o <C -^-j x j^- Under this condition, the 
relation between the parameters of the diffusion model and the equivalent leaky model follows from 
the equation X^f = X^ that leads to the expression for diffusion coefficient 



D = ^cH/(2X lb ), (3-8) 

where (3 = v/c. In this simplified model the path length distribution is a decreasing exponential 
function with a mean path length dependent on energy. Such a dependence is shown in equation 
(3-9)) where the parameters Xq, Rq, and a are determined by the density of the matter disk, the 
size of the halo and the spectrum of magnetic turbulence that sets the diffusion coefficient's rigidity 
dependence. 

X w = X p g/cm 2 at R < i? GV, X w = X (3(R/ R GV)~ a g/cm 2 at R > R GY (3-9) 

where R is the particle rigidity and (3 = v/c. This parameterization is valid for particles in the 
interstellar medium with energies from about 0.4 GeV/n to 300 GeV/n where data on secondary 
nuclei are available. It is worth noting that the same escape length (equation (3-9)) satisfactorily 
reproduces both B/C and Sub-Fe/Fe ratios (no need for path length "truncation"). It should be 
borne in mind, however, that this diffusion model does not predict the sudden change of the mean 
grammage as a function of energy that is displayed in equation (3-9) at a rigidity of Rq GV or any 
other rigidity. This is a strictly (though widely accepted) ad hoc construction that appears to be 
required by the data. We will see that the models following this one do not require this behavior 
to be externally imposed. 

The physical interpretation of the empirical equation (3-9) for R > Rq GV can be given in the 
framework of the diffusion model by referring to equation (3-8). The theory of particle resonant 
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scattering and diffusion in the turbulent interstellar medium predicts the scaling of the diffusion 
coefficient D res = n(3R a , where constant k is determined by the level of hydromagnetic turbulence 
with the spectrum Wkdk oc k~ 2+a dk, a =const. The scattering of particles with Larmor radius r g is 
mainly due to the interaction with inhomogeneities of the scale 1/k ~ r g . Since the time to diffuse 
to the galactic boundary, and hence the mean grammage traversed, is inversely proportional to 
the diffusion coefficient the interpretation is clear. The observations of interstellar turbulence are 
consistent with the existence of a single power-law spectrum with 0.2 < a < 0.6 at wave numbers 
10 _20 cm _1 < k < 10 _8 cm _1 , see Ruzmaikin et al. (1988). The Kolmogorov spectrum, usually 
considered to be representative of the interstellar spectrum corresponds to a = 1/3. 



4. Galactic wind model. 

We consider a one-dimensional galaxy model shown in Figure (2) (the coordinate z is per- 
pendicular to the galactic plane). Cosmic ray sources and interstellar gas are concentrated in a 
thin disk at z = (For details of this model see Jokipii (1976); Jones (1979)). The distribution 
function f(z,p) obeys the equation 



+ u— + —5{z)f + —— t 

dz dz dz dz?> dp m p 2 dp 

where u is the wind convection velocity. 



q (p)5(z). (4-1) 
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Fig. 2. — Simplified galactic model with a halo wind; galactic matter disk is of infinitesimal thick- 
ness. 
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We shall assume that the wind velocity u is constant and directed outward from the galactic 
plane. There is a cosmic ray halo boundary at i z = H where cosmic rays freely exit from the 
Galaxy Using the same prodecure as in section (3) we integrate equation (4-1) in the vicinity of 
galactic plane (lim e ^o fl e dz(...)) one can find the boundary condition at z = + e : 



ori <9/b Zupdfo fiva 1 d 

~ 2V ~x ^T~x~ 1° + ~2 ^~ 

oz 6 op m p z op 



' 2 /A) , ' 
P fo 

m 



= Qo(p), 



(4-2) 



where fo(p) = f(z = 0,p) is the distribution function at the galactic midplane. 
The solution of equation (4-1) at z ^ under the boundary condition /(i z \= H) = is: 

1 -exp(— u(H- i z \)/D) 



f = h : 



1 - exp (-uH/D) 



(4-3) 



Calculating —D^fr from equation (4-3) and substituting it in equation (4-2), one can get the 



following equation for fo(p) '■ 
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Equation (4-4) gives the following equation for the cosmic-ray intensity I{Ej t ) 

'dE k 
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dE k 



dx J ad 



dx 



2_j = qq(.p)p 2a 

m fiv 



with the effective escape length 



X w = g (1 - exp (-uH/D)) , 



and the adiabatic energy loss rate per g/cm 2 

fdE k \ 



2u 



\ dx ) ad 3 ^ C 



- — \J E k (E k + 2E Q ). 



(4-5) 



(4-6) 



(4-7) 



At this point we have arrived at an equation that is essentially the same as for the thin disk-halo 
model in the sense that the path length distribution will be an exponential with the mean path 
length given by X w . We may therefore apply the modified weighted slab technique to solve the 
equations. 



As before, the asymptotic scaling of X w in this case is 

LLV 

X w = — oc v at small rigidities (when uH/D » 1), 

LI 



(4-8) 
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and 

X w = oc R~ a at large rigidities (when uH/D < 1).. (4-9) 

It should be noted that here the flattening of the curve of mean grammage vs. energy at low energy 
arises naturally from the model and does not need to be inserted by hand. 

Equation (4-6) may be presented in the following form useful to fit the observations: 

X.=0X c (l- ap (- 7 ^ jt )). (4-10) 

Here Xq = (fic)/(2u), and Rq = (uH/ kq) 1 ^ 11 . The adiabatic energy loss term, equation (4-7) in 
these notations reads as 

It is worth noting that here we consider only a simple model with constant wind velocity. In 
more realistic models the wind velocity depends on th distance from the galactic plane, e.g. in a 
self consistent model of a cosmic ray driven galactic wind (Ptuskin et al. 1997). 



5. Turbulent diffusion 

In this model there is no regular convective (wind) transport, rather the flow is turbulent and 
the transport of the particles is more random. In fact the convection in this flow is of the same 
nature as diffusion with the diffusive properties determined by the random flow. (See Figure (3)) 
Equation (4-1) becomes 



d r^df Liva . , Id 
oz az m p z op 



P 



dp 
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q (p)S{z). 



(5-1) 



The equation for I{Ek) is the following 
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(5-3) 
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Fig. 3. — Turbulent diffusion model, halo is diffusive with two types of diffusion. 

. We assume now that cosmic ray diffusion is provided simultaneously by turbulent diffusion 
with the diffusion coefficient Dt that does not depend on particle energy (may be estimated as 
D t = u t L t /3, where ut and L t are the characteristic random velocity and correlation scale of 
large-scale turbulent motions of the interstellar gas) and by resonant diffusion with the diffusion 
coefficient D res = [3KoR a (here Ko = const, a = const) provided by the scattering on hydromagnetic 
turbulence as was discussed in the previous section. The total diffusion coefficient which appears 
in equation (5-3) is equal to 

D = D t + D res . (5-4) 



Equation (5-3) may be presented in the following form that is useful to fit the observations: 

Xdlf = 1 + (3{R/R f (5 ~ 5) 



Here X = (ficH)/(2D t ), and R = (A/«o 



We note again the flattening of the mean grammage vs. energy curve for low rigidities, as 
in the previous wind model, arises naturally from the physical picture. Actually the reason is 
the same in both models; below a particular rigidity kinetic diffusion becomes slower than the 
convective transport in removing particles from the galaxy and the particles residence time becomes 
independent of rigidity or energy. 
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6. Stochastic Reacceleration. 



This model (Seo & Ptuskin 1994) assumes no convective or wind motion in the system( see 
Figure (4)). As before, the spatial diffusion is provided by scattering on random hydromagnetic 
waves with a spectrum that results into D oc vR a . For a Kolmogorov spectrum of turbulence 
a = 1/3; but we shall consider the value of a as a free parameter, in the same way as in the models 
discussed above. This determines the behavior of the escape length X e = XoR~ a (see equation (4- 
9)) at all energies. So far this does not differ from the diffusion model but here we consider the fact 
that the turbulence is not static, rather, the magnetic fluctuations move with the Alfven velocity 
producing a diffusion in momentum as well as in space. Stochastic acceleration with a diffusion 
coefficient in momentum K ~ p 2 v\/D (va is the Alfven velocity) essentially modifies the spectra of 
primaries and secondaries below about 10 GeV/n and can produce the characteristic peak in B/C 
ratio at few GV. The acceleration becomes inefficient at high energies and the model is reduced 
to a simple diffusion model without reacceleration and with the escape length X e = XoR~ a at 
E > 20 - 30 GeV/n. 
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Fig. 4.— 



Reacceleration model, halo is diffusive in momentum as well as space. 
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We use the same notations as Seo and Ptuskin and use their equation (12) for the calculations: 
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Here the parameter a, which is defined as 



32 



a 



h a ( v a 



(6-1) 



(6-2) 



3a(4 - a 2 )(4 - a) H \fic / 

(its dimension is (g/cm 2 ) -2 ), determines the efficiency of reacceleration, h a is the height of reac- 
celeration region. In this work we have taken h a /H = 1/3. 

The parameters we have to find from fitting the data are Xq, a and a since we will not prescribe 
the spectrum of the interstellar turbulence. 



7. Fits to data 

We fit the four models discussed above to a collection of data compiled by Stephens & Streit- 
matter (1998) for the B/C and for the Sub-Fe/Fe ratios. We determined which parameters best 
fit (in the weighted least squares sense) both ratios simultaneously for each model. The parameters 
found for each model are given in Table (1). It should be noted that because considerable comput- 
ing is required for each set of parameter values the search in parameter space was not automated. 
It was performed by hand and thus we can not guarantee that the fits that we have found are 
rigorously least squares fits, they are simply the smallest values we could find in our searches. 



Fitted Parameters 



Model 


X (g/cm 2 ) 


Ro (GV) 


a 


X 2 (normalised) 


Disk-Halo Diffusion 


11.8 


4.9 


0.54 


1.3 


Wind 


12.5 


11.8 


0.74 


1.5 


Turbulent Diffusion 


14.5 


15.0 


0.85 


1.8 


Stochastic Reacceleration 


9.4 


a = 2.6 x 1(T 3 


0.30 


1.8 



Table 1: Best fit parameters for fitting secondary /primary ratios for the models considered. For 
Minimal Reacceleration there is no Rq parameter, rather the strength of acceleration is given by 
the dimensionless parameter a. 



These parameters imply physical quantities for the different models given in Table (2). Notice 
that certain parameters are meaningful for specific models only, not all models in general. 
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Physical Parameters 



Model 


D res (cm 2 /s) 


D t (cm 2 /s) 


u (km/s) 


V a (km/s) 


Disk-Halo Diffusion 


2.0 x 10 28 /? 








Wind 


7.2 x 10 27 /? 




29 




Turbulent Diffusion 


3.8 x 10 27 /? 


3.8 x 10 28 






Stochastic Reacceleration 


5.9 x 10 28 /3 






40 



Table 2: Physical parameters implied by the fit parameters of Table 1. The values of D res should 
be multiplied by R a with R in GV and a taken from Table 1 



In the turbulent diffusion model the turbulent diffusion coefficient can be estimated as Dt = 
u t L/3, where ut is the characteristic velocity and L is the characteristic scale of the turbulent 
motions. Since we would not expect the magnitude of the velocity ut to exceed 100 km/s that leads 
to the very large value of L > 0.76H. Thus this model of cosmic ray transport requires a value of 
D t that is difficult to reconcile with acceptable values of parameters of the interstellar turbulence. 

These fits are displayed with the data in Figures (5) and (6). After determining the best fit 
parameters we then used them to propagate primary spectra for C, O and Fe as discussed in Section 
(8). If one compares Figures (5) and (6) with Figures (7) and (8) in Section (8) one can see how 
much more sensitive secondary to primary ratios are to the propagation model than are simple 
primary spectra. Although the models were all chosen to fit the same secondary to primary data 
one can still see significant differences in the model curves whereas in the plots of carbon and iron 
the model curves are quite hard to distinguish from one another. 

8. Source spectrum. 

Using the above determined propagation models, we tested different source spectra to fit 
observed spectra of primary nuclei C, Fe at energies 0.5-100 GeV/n. The best fit for the disk-halo 
diffusion, wind and turbulent diffusion model is provided by the source spectrum Q oc i?~ 2 35 ; for 
the minimal reacceleration model the best fit was obtained with Q oc (3R~ 2A0 Figures (7) and (8) 
show how the different propagation models fit the observations assuming these source spectra. One 
can expect that asymptotically at very high energies, E > 100 GeV/n, the relation 7 = 7 S + a, 
where 7(73) is the exponent of the observed (source) differential spectrum 1(E) oc E^" 1 ) is fulfilled. 
It is of interest that three of the four models were well fit with a source spectrum that was a 
simple power law in rigidity. Only the minimal reacceleration model required a modification at low 
(non-relativistic) energies. A similar behavior was found by Heinbach & Simon (1995) in their 
investigation of a reacceleration model. 
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Fig. 5. — The least squares fit to observed B/C ratios, in four propagation models: turbulent 
diffusion (dashed lines), wind (dotted lines), reacceleration (dash-dotted lines), and the disk-halo 
diffusion model with the diffusion coefficient given by equation (3-8) (solid lines). is the force 
field approximation solar modulation parameter. Data are from a comprehensive compilation by 
Stephens & Streitmatter (1998). 

9. Anisotropy constraint 

The observed galactic cosmic rays are highly isotropic. The amplitude of the first angular 
harmonic of the cosmic ray distribution in the interstellar medium near the solar system is ap- 
proximately equal to 5 = (0.5 — 1.0) x 10~ 3 at energies 10 12 — 10 14 eV (Nagashima et al. 1989; 
Cutler & Groom, 19991; Aglietta et al. 1993; Alekseenko et al. 1993). The anisotropy is almost 
energy-independent but a slow increase with energy or even irregular behavior at energies close to 
10 14 eV are not excluded. 

Possible interpretation of these observations is based on a contribution of local cosmic ray 
sources (Dorman et al. 1984). Supernovae and their remnants are assumed to be the instantaneous 
point sources. The cosmic ray density is determined by the total contribution from numerous 
sources (SN outbursts occur in about each 30 yr, whereas the cosmic ray confinement time in the 
Galaxy exceeds 10 7 yr). However, the anisotropy may be defined by an individual nearby source. 
If the diffusion coefficient D(E) increases with energy, the contribution of an individual outburst 
at distance r gives rise to an anisotropy which depends nonmonotonicly on energy. A source's 
contribution to the anisotropy amplitude reaches its maximum when t ~ r 2 /D, where t is the 
age of instantaneous source. Analysis of the list of supernova remnants and pulsars indicates that 
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Fig. 6. — The least squares fit to observed Sub-Fe/Fe ratios, in four propagation models: turbulent 
diffusion (dashed lines), wind (dotted lines), reacceleration (dash-dotted lines), and the disk-halo 
diffusion model with the diffusion coefficient given by equation (3-8) (solid lines). is the force 
field approximation solar modulation parameter. Data are from a comprehensive compilation by 
Stephens & Streitmatter (1998). 



Geminga, Vela, Lupus Loop, Loop III, and some others may prove to be the sources sustaining the 
anisotropy observed at 10 12 -10 14 eV. 

One has to check the compatibility of the propagation models discussed in the previous sections 
with the data on cosmic ray anisotropy. Each of these models assumes a specific energy depen- 
dence of cosmic ray leakage from the Galaxy. It might result in an anisotropy which exceeds the 
observational limit. Effect of the solar wind masks the anisotropy of galactic comic rays at energies 
less than about 10 12 eV for an observer at the Earth. Below we assume that the dependence of 
cosmic ray diffusion on energy determined from the observations of secondary to primary ratios in 
cosmic rays at energies up to about 10 11 eV can be extrapolated to energies 10 12 — 10 14 eV. This 
simplifying assumption is based on the observation that there is no any drastic change in the total 
cosmic ray energy spectrum in the energy range 10 10 - 10 14 eV that would be indicative of change 
in energy dependence of cosmic ray transport. 

The equation for the amplitude of cosmic ray anisotropy perpendicular to galactic plane is the 
following (see e.g. Berezinskii et al. (1990)): 



A 3 fn d f m P d f\ 
S = -^f{ D d-z +U 3^) 



(9-1) 
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Fig. 7. — The least squares fit to the C spectrum in the same propagation models as indicated in 
Figures (5) and (6). 

which includes both the diffusion and convection fluxes of cosmic rays. 

It is easy to show that diffusion dominates over convection at high enough energies E > 100 
GeV and that effects of ionization energy losses and reacceleration are not essential at these energies 
in the models we discuss here. Notice also that cosmic ray anisotropy is mainly determined by the 
most abundant proton component of cosmic rays that is subject to insignificant effect of nuclear 
interaction with the interstellar gas. In these conditions, the second term in brackets in the last 
equation can be omitted, and the cosmic ray transport equation can be presented in a simple form 
as 



Now it is easy to show that the anisotropy for an observer just above the <5-plane with cosmic ray 
sources is approximately equal to 



in all diffusion models . Here X is the escape length which obeys the equation X m nflcH/ (2D) at 
high particle energies. 

Cosmic ray anisotropy inside the source region is smaller than Sq. For cosmic ray sources 
uniformly distributed through the disk with a total thickness 2h, the anisotropy at distance z 
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Fig. 8. — The least squares fit to the Fe spectrum in the same propagation models as indicated in 
Figures (5) and (6). 

from the central galactic plane is estimated as S z ~ S^z/h at \z\ < h (Ptuskin 1997) . With 
parameters found in the previous sections and extrapolated to energy 10 14 eV and with assumption 
that z/h = 0.1 (z = 20 pc, h = 200 pc), we have the following expected values of the anisotropy 
5 2 (10 14 eV) which is due to the streaming of cosmic rays perpendicular to galactic plane: 7 x 10~ 3 
in the basic diffusion model, 4 x 10~ 2 in the model with turbulent diffusion, 2 x 10~ 2 in the galactic 
wind model, 1 x 10~ 3 in the minimal reacceleration model. These values are uncomfortably large, 
even for the minimal reacceleration model 

Of course, the models of the galaxy that we have employed here are highly simplified, primarily 
in the high degree of symmetry and smoothness that they exhibit. This means that the anisotropics 
that we have calculated are most likely lower limits to those that would be obtained from more 
complex models. It is always possible to construct models with the solar system in a rather favored 
position of symmetry that would produce a smaller value of the anisotropy but such models, lacking 
any particular knowledge of their reality, are a priori unlikely. 



10. Conclusion. 

The objective of the present work was to investigate several models of cosmic ray propagation 
and nuclear fragmentation in the interstellar medium using the most recent set of spallation cross 
sections and employing the modified weighted slab method which allows us to obtain exact solutions. 
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It is clear from a comparison of Figures (5) and (6) on the one hand and Figures (7) and (8) on 
the other that the secondary to primary ratio is a much more sensitive function of the propagation 
model (and relatively insensitive to the assumed source spectrum) than is the relation of an observed 
spectrum to the source spectrum. It is for just this reason that we used this ratio to test our various 
models before then deducing primary spectra. The standard flat-halo diffusion model, the models 
with turbulent diffusion, the model with constant wind velocity, and the diffusion model with 
reacceleration were tested. For each model we then picked the parameter set that best fit the data. 
The primary spectra were then calculated for each model using these parameters. All these models 
are able to explain the decrease of secondary/primary ratios at energies below a few GeV/n. The 
turbulent diffusion and the wind transport work simultaneously with resonant diffusion. The last 
process dominates in cosmic ray transport at high energies. In order to reproduce a sufficiently sharp 
bend in the secondary /primary ratio in these models the diffusion must have a strong dependence 
on rigidity D res oc f3R a , where a = 0.74 - 0.85, at least at low energy (E < 30 GeV/n). However, 
this could lead to a problem. The extrapolation of such strong rigidity dependence of diffusion 
to energies 10 3 — 10 5 GeV produces anisotropies that are at strong variance with the observed 
anisotropy at these energies. 

The disk-halo diffusion model has less severe problem with the explanation of observed small 
anisotropy than the models with turbulent diffusion and wind. Also, the scaling of diffusion co- 
efficient on rigidity D res oc (3R - 54 at large rigidities R > 4.5 GV in this model is probably not 
in contradiction with the observations of interstellar turbulence. At the same time, the sharp 
change-over of the diffusion to the regime D = const at rigidities R < 4.5 GV can not be natu- 
rally explained by the kinetic theory of particle scattering in random magnetic fields. The ad hoc 
elucidation of the observed peak in the secondary-to-pramary ratios makes the model implausible. 

The reacceleration model is favored as regards to its natural description of the energy depen- 
dence of particle transport in the Galaxy with a close to Kolmogorov spectrum of turbulence since 
the scaling D res oc (3R - 3 at all rigidities is assumed in this case. The structure of the peak in the 
secondary-to-primary ratios is reproduced in this model. The model with reacceleration has no 
problems with the account for low anisotropy of cosmic rays. However, the predicted weak energy 
dependence of secondary-to-primary ratios at energies > 20 GeV/n is not in agreement with the 
high energy HEAO-3 data by Binns et al. (1981) on Sub-Fe/Fe ratio, see Figure 6. 

The source spectra for primary nuclei derived in all propagation models are close to R~ 2 - 3 — 
R~ 2A that is uncomfortably steep compared with the predictions of supernova shock acceleration 
models (Baring et al. 1999; Berezhko &; Ellison 1999). Strictly speaking, this refers to the limited 
energy range up to about 30 GeV/n where statistically accurate data on secondary and primary 
nuclei are available. Even if one does not assume that a single power law source spectrum continues 
to much higher energies, the soft depends of diffusion on energy in the reacceleration model imply 
a steep source spectrum close to R~ 2A to fit the observed spectrum of primaries close to i?~ 2 ' 75 . 
This problem may be not so pressing for the disk-halo diffusion model. 
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It is clear that measurements of secondary to primary ratios at higher energy are needed since 
it is here that the models begin to seriously diverge. Such future missions such as ACCESS would 
be very helpful in determining which pictures of galactic cosmic-ray propagation are feasible and 
which are not. 
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Figure Captions 

Fig. 1. — Simplified galactic model with a diffusive halo, galactic matter disk is of infinitessimal 
thickness. 

Fig. 2. — Simplified galactic model with a halo wind; galactic matter disk is of infinitessimal 
thickness. 

Fig. 3. — Turbulent diffusion model, halo is diffusive with two types of diffusion. 

Fig. 4. — Reacceleration model, halo is diffusive in momentum as well as space. 

Fig. 5. — The least squares fit to observed B/C ratios, in four propagation models: turbulent 
diffusion (dashed lines), wind (dotted lines), reacceleration (dash-dotted lines), and the disk-halo 
diffusion model with the diffusion coefficient given by equation (3-8) (solid lines). is the force 
field approximation solar modulation parameter. Data are from a comprehensive compilation by 
Stephens & Streitmatter (1998). 

Fig. 6. — The least squares fit to observed Sub-Fe/Fe ratios, in four propagation models: turbulent 
diffusion (dashed lines), wind (dotted lines), reacceleration (dash-dotted lines), and the disk-halo 
diffusion model with the diffusion coefficient given by equation (3-8) (solid lines). is the force 
field approximation solar modulation parameter. Data are from a comprehensive compilation by 
Stephens & Streitmatter (1998) 

Fig. 7. — The least squares fit to the C spectrum in the same propagation models as indicated in 
Figures (5) and (6). 

Fig. 8. — The least squares fit to the Fe spectrum in the same propagation models as indicated in 
Figures (5) and (6). 



